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RATIONALE 

Phytoplankton  directly  influence  many  optical  and  acoustic 
characteristics  of  ths  oceanic  water  column.  Under  nutrient  sufficient 
conditions  (Falkowsd  et  al.,  1992),  phytoplankton  activity  is  regulated 
both  by  the  varyin*  photosynthetically  active  photon  fluence  density 
(PPFD)  exposure  and  by  the  photosynthetic  potential  of  the  individual  cells 
throughout  the  wat<  r  column  (Maira,  1978a;  1978b). 

PPFD  at  the  sea  surface  varies  over  a  range  of  time  scales  due  to  clouds, 
and  day/night  and  seasonal  cycles  (Kirk,  1983).  Subsurface  PPFD  exhibits’ 
spectral  changes  w  th  depth  due  to  wavelength-dependent  exponential 
attenuation  (Kirk,  1  ?83),  Additional  temporal  variation  in  subsurface  PPFD 
results  from  dissolved  organic  matter  (DOM)  and  particulate  organic  matter 
(POM)  patchiness  (Parsons  et  al.,  1984a).  Near  surface  water  motion,  the 
net  effect  of  turbulmt  mixing,  surface  and  internal  gravity  waves,  wind 
drift,  Langmuir  circ  ilation  and  larger  scale  flow  (Denman  and  Gargett, 
1983),  transports  al|al  cells  through  the  fluctuating  PPFD  field.  Some  cells 
act  like  neutral  La^rangian  particles,  but  most  possess  an  inherent  motility 
due  to  buoyancy  (Snayda,  1970;  Walsby  and  Reynolds,  1980)  or 
swimming  (Roberts,  1981;  Kamykowski  et  al.  1988).  The  physical  and 
biological  motion  vectors  combine  to  determine  phytoplankton  trajectories. 
Since  the  emphasis  here  is  PPFD  exposure  of  individual  phytoplankton 
cells,  vertical  displacement  is  of  primary  concern  (Kamykowski,  1990; 
Yamazaki  and  Kamykowski,  1990).  One  aspect  of  the  work  in  his  report 
dealt  with  the  deve  opment  of  the  Physical  Motion  Instrument  (PMI)  to 
monitor  water  motitin  over  biologically  significant  time  scales. 

In  a  previous  review,  Falfcowski  (1984)  emphasized  the  need  for  in -situ 
work  that  includes  both  O2  and  CO2  fluxes  to  advance  our  understanding  of 
the  time-space  relationships  that  influence  phytoplankton  photosynthesis. 
An  in-situ  approach  includes  the  complex  relationships  between  PPFD 
intensity  and  the  li|;hi  spectrum  in  the  natural  water  column.  The 
combined  measurement  of  O2  and  CO2  fluxes  provides  a  useful  indicator  of 
the  instantaneous  stite  of  photosynthesis  in  response  to  fluctuating  light 
fields.  Although  wcrk  with  diverse  natural  communities  provides  useful 
information,  an  acc  irate  assessment  of  how  distinct  phytoplankton  species 
respond  to  fluctuati  ig  light  benefits  from  work  with  axenic,  single-species 
cultures  under  realistic  upper  ocean  conditions.  Various  aspects  of  a 
species’  complex  pi  otosynthetic  mechanism  (Falkowski,  1984)  result  in 
both  rapid  response;  due  to  minor  biochemical  changes  accompanying 
fluorescence  modifii  ation  and  lengthier  responses  such  as  sun-shade 
adaptation  (Prezelin,  1981)  and  photoinhibition  (Neale,  1987).  When  the 
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environmental  variability  exceeds  the  adaptability  of  a  given  species, 
succession  to  a  nev'  species  with  a  more  appropriate  range  of 
photosynihetic  poteitial  often  occurs  (Sakshaug  et  al.,  1987;  Caron  et  al., 
1988).  Therefore,  measurements  on  single  species  from  representative 
functional  groups  (I ‘ichardson  et  al.,  1983)  must  be  analyzed  to  examine 
the  scope  of  intersi  ecies  responses.  A  second  aspect  of  this  report  deals 
with  improvements  to  the  ONR-funded  Self-contained  Underwater 
Photosynthesis  App  iratus  (SUPA),  a  facility  that  originally  measured  the 
temperature,  PPFD  exposure,  and  the  O2  and  CO2  fluxes  of  an  enclosed 
phytoplankton  cultire  at  sea 

To  sort  out  tl  e  relative  importance  of  various  processes  theoretically, 
scaling  procedures  ypically  are  utilized.  Lewis  et  al.  (1984a;  1984b)  place 
the  combined  effect  of  varying  PPFD  and  photosynthetic  potential  of 
existing  phytoplankton  populations  at  a  given  time  into  a  robust 
nondimen  si  onal  framework  based  on  simple  assumptions  concerning 
turbulent  mixing,  eiphotic  zone  depth  and  photoadaptation  rate.  To  predict 
future  states  of  productivity  over  regions  of  the  ocean  from  current 
phytoplankton  distributions,  models  are  require  with  appropriate 
parameterizations  oj  important  small  scale  physical  and  biological 
processes.  Both  EuUrian  (Martin,  1986;  Price  et  al.,  1986)  and  Lagrangian 
(Woods  and  Onken,  1983;  Kamykowski,  1990;  Yamazaki  and  Kamykowski, 
1990)  models  of  particle  dispersion  in  the  upper  ocean  are  useful  in 
different  contexts  ( -ande  and  Lewis,  1989).  Continued  real  progress, 
however,  requires  field  data  to  test  model  assumptions  and  provide  insight 
into  how  dynamics  of  relevant  vertical  water  motion  (Weller  and  Price, 
1988;  Gargett,  1985)  and  organism  behavior  affect  the  kinetics  of 
photoadapiation  (Cullen  and  Lewis,  1988;  Kamykowski,  1979;  Denman  and 
Marra  et  al.,  1985;  Neale  and  Richardson,  1987)  and  the  photosynihetic 
potential  of  individual  species  (Richardson  et  al.,  1983).  Such  field 
observations  require  innovative  instrumentation  like  PMI  and  SUPA  to 
determine  the  cohe  rences  among  the  photosynthetic  response,  vertical 
water  motion  and  incident  radiation.  A  third  aspect  of  this  report  is  to 
decribe  models  thai  combine  water  motion,  behavior  and  physiology  using 
computer  facilities  ihat  were  expanded  using  the  ONR-SORP  funds.  These 
models  were  generated  with  the  expectation  that  supporting  data  from 
PMI  and  SUPA  eventually  will  be  available. 


GENERAL 

This  ONR-SORP  grant  was  written  to  support  collaborative  work  with 
Dr.  Tom  Curtin,  ONR.  Separate  ONR  and  NSF  funding  was  available  to 
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develop  and  deploy  the  Self-contained  Underwater  Photosynthesis 
Apparatus  (SUP A)  discussed  below  and  to  model.  The  funds  in  this  ONR- 
SORP  grant  were  piimarily  allocated  to  support  the  development  of  the 
Physical  Motion  Instrument  (PMI)  by  Dr.  Curtin  to  provide  information  on 
how  SUPA  should  b  ?  moved  in  the  upper  mixed  layer  of  the  ocean.  Some 
funds  also  were  allocated  to  upgrade  SUPA  and  to  the  expedite  the 
modeling  effort. 


PHYSICAL  MOTION  INSTRUMENT  (PMI) 

The  PMI  is  an  array  of  current  sensors  nested  symmetrically  along  two 
orthogonal  axes  each  one  meter  long  erossing  at  the  eenter.  When  the 
normal  to  the  plane  is  oriented  vertically,  the  vertical  velocity  component 
and  one  horizontal  component  are  measured  independently  at  the  end  of 
each  axis.  At  the  center,  the  two  horizontal  velocity  components  are 
measured  at  0.4  cm  intervals  over  a  range  of  25  cm  along  each  axis.  This 
arrangement  permits  calculation  of  velocity  gradients  and  thus  the  vertical 
component  of  vorticiiy  over  a  range  of  scales,  and  enables  comparison  of 
the  small  scale  struc  ture  with  the  one  meter  scale  flow  gradients  derived 
from  the  comers.  Findamental  issues  are  the  spatial  and  temporal 
coherence  of  the  ve!  ocity/vorticity  field  in  this  scale  range  and  the  best 
sampling  strategy  fer  measuring  ocean  velocity  fluctuations  that  affect 
phytoplankton  productivity. 

Performance  criteria  considered  in  evaluating  the  best  flow  sensing 
method  for  the  objec  tives  included:  (1)  low  zero  velocity  drift,  such  that  a 
mean  flow  can  be  accurately  measured  in  both  high  amplitude  unsteady 
flow  and  low  amplitude  quiescent  conditions;  (2)  adequate  temporal 
response  to  resolve  high  frequency  surface  gravity  wave  motion;  (3)  linear 
amplitude  response  ever  the  full  range  of  velocity,  typically  0  to  2  m/s;  (4) 
cosine  or  otherwise  stable  angular  response  in  the  plane  of  measurement; 
(5)  negligible  resporse  to  flow  orthogonal  to  the  plane  of  measurement 
(cross  talk);  (6)  low  sensitivity  to  fouling  to  minimize  calibration 
compromise  over  lo  ig  intervals  of  deployment;  (7)  mechanical  ruggedness 
to  withstand  deploy]  nent  shocks,  storm  event  forces  and  long  term  wear; 
(8)  potential  for  coincident  temperature,  conductivity,  pressure  and 
particle  type/concen  ration  measurement;  (9)  low  power  consumption  to 
enable  long  term,  remote  deployment;  (10)  low  cost  and  modular  design  to 
enable  large  arrays.  Non-mechanical  speed  and  direction  sensing 
components  best  saiisfy  frequency  response  and  ruggedness  criteria.  The 
principal  types  of  i  on-mechanical  speed  sensors  are  electromagnetic, 
acoustic  and  laser.  Electromagnetic  sensors  are  based  on  the  emf  induced 
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by  a  conductor  (seavater)  moving  through  a  magnetic  field  (earth's  or 
locally  produced).  Acoustic  and  laser  sensors  are  based  either  on 
backscatier  doppler  shift  or  dual  pulse  correlation,  forward  scatter  dual 
pulse  coherence  or  leciprocal  time  of  flight.  The  electromagnetic  (EM) 
method  was  chosen  to  minimize  the  sensitivity  of  the  measurement  to 
suspended  particle  ,ind  bubble  concentrations  and  non-Lagrangian 
behavior. 

EM  sensors  in  various  compact  configurations  have  been  fabricated 
previously,  and  many  have  worked  acceptably.  A  central  problem  with 
this  type  of  sensor  !ias  been  the  low  induced  output  voltage  when  power  is 
limited  for  magnetic  field  generation.  The  low  signal-to-noise  ratio 
necessitates  high  gain  amplification  leading  to  problems  with  zero  drift.  In 
addition,  the  emf  simsing  electrodes  have  inherent,  transient 
electrochemical  poteitial  differences  as  well  as  potential  differences 
related  to  local  bou  idary  layer  modification.  The  relatively  slow  drifting 
electrochemical  differences  can  be  effectively  eliminated  by  ac  driving, 
either  electrically  oi  mechanically.  The  boundary  layer  related  differences 
can  be  controlled  b}  physically  isolating  electrodes  from  the  flow  and/or  a 
geometry  with  reprcducable  boundary  layer  effects  related  to  the  flow 
field  measurement  volume  (the  e-folding  volume  of  magnetic  field 
intensity).  Examples  and  discussion  of  boundary  layer  effects  have  been 
given  by  Aubrey  ei  al.  (1984).  Complete  calibration  includes  potential 
differences  due  to  loth  flow  through  the  measurement  volume  and  flow 
altering  the  electrote  boundary  layer. 

The  PMI  incoi]  orates  two  types  of  EM  sensor:  macroscale  velocity 
sensors  (MVS)  at  thj  ends  of  each  axis  and  small-scale  velocity  sensors 
(SVS)  along  the  axe: .  For  the  MVS,  a  new  low  cost  amplifier  circuit  was 
designed  for  low  si;  ;nal-to-noise  ratio  and  high  common  mode  rejection, 
providing  very  low  tero  drift,  and  a  new  geometry  was  developed  to 
control  boundary  la  ^er  and  fouling  problems  and  improve  ruggedness  and 
deployment  flexibility.  The  magnetic  field  is  generated  by  a  Helmholtz  coil 
arrangement  as  suggested  by  Olson  (1972),  however  the  electrodes  and 
coils  are  constructed  differently.  The  coil  units  are  solid,  rather  than  open, 
providing  a  ducted  flow  measurement  volume  within  a  highly  uniform 
field.  Identical  sets  )f  electrodes  are  flush  mounted  on  opposing  faces  of 
the  sensor  (co-plane/  with  the  coils)  and  connected  in  parallel.  The 
complete  configuration  is  mechanically  robust  and  hydrodynamically 
symmetric,  providin;'  reliable,  precise  measurements  (Sorrell,  ei.al.,  1990). 

The  MVS  has  been  designed  to  incorporate  measurement  of  tilt, 
temperature,  conductivity,  pressure  and  optical  beam  attenuationwithout 
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extern^  modificEtion  to  its  goometry.  For  vajisblos  such  es  conductivity 
and  bcEm  Ettenuatioi  the  configuration  is  ideal.  Another  design  feature  is 
the  sensor's  function  as  an  integral,  load  bearing  element  within  a  semi¬ 
rigid  spar.  A  modulsT  system  may  be  configured  as  a  series  of  discrete 
sensor  units  interco  mected  mechanically  by  variable  length  spar  elements 
and  electrically  by  in  interface  bus  to  a  central  compass/controller/data 
logger/power  supply.  The  result  is  a  compact,  integrated  sensor  array  that 
can  be  networked  for  cost-effective,  long  term  measurement  of  gradients. 

The  SVS  consists  of  two  orthogonal  linear  arrays,  each  with  64 
electrodes  (0.4  cm  uniform  spacing),  flush  mounted  on  an  open  face  inside 
and  co-planer  with  a  single  coil  creating  a  dipole  magnetic  field.  The  flow 
induced  electrical  potential  difference  between  any  two  pairs  of  electrodes 
may  be  calculated  f  om  the  measured  emf  across  the  arrays,  enabling  the 
coherence  over  a  range  of  scales  to  be  determined  from  time  series  data. 
Data  obtained  during  field  trials  of  the  PMI  at  Duck,  NC,  have  proven  the 
ability  of  the  SVS  to  detect  signals  in  the  ocean  even  under  benign,  low 
energy  conditions.  Further  laboratory  work  is  necessary,  however,  to 
precisely  define  the  transfer  function  between  measured  voltage  gradients 
and  the  structure  of  the  inducing  flow  field. 


SUPA 


A  Self-containt'd  Underwater  Photosynthesis  Apparatus  (SUPA)  has 
been  designed  to  mr.asure  photosynthetic  response  of  a  single  species 
phytoplankton  cultun;  to  PPFD  variations  in-situ  (Kirkpatrick,  1990).  The 
PMI  has  been  desigied  to  measure  velocity  field  temporal  fluctuations 
over  a  range  of  spat  al  scales  simultaneously  from  1  m  to  0.4  cm  to 
examine  the  coherence  of  velocity  and  velocity  gradient  (e.g.,  vorticity)  in  a 
local  volume.  Couple!  together  with  SUPA,  PMI  can  also  be  used  as  a 
vertical  motion  sense  r  for  interactive  control  of  SUPA  operating  in  a  quasi- 
Lagrangian  mode.  The  prototype  instruments  have  been  successfully 
demonstrated  in  field  trials  (Kirkpatrick,  et  al.,  1990)  and  currently  are 
being  refined. 

SUPA  resolves  photosynthetic  and  physiological  responses  to  in  situ 
PPFD  and  temperatire  variations  by  simultaneous  measurement  of  both 
carbon  uptake  and  cxygen  evolution  at  one  minute  intervals.  Carbon 
assimilation  is  deter  nined  through  the  measurement  of  pH  (Siumm  and 
Morgan,  1981;  Axehson,  1988),  alkalinity  corrections  due  to  uptake  of 
nitrate  (Goldman  an  1  Brewer,  1980).  Measurements  include:  the  minutes 
to  hours  response  ol  a  culture  to  light  fluctuations,  the  PPFD-dependent, 
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light-saturated  photosynthetic  rates  and  inhibition  parameter,  the 
variation  of  dark  respiration  rate,  and  the  PPFD-dependent  variation  in  PQ, 
an  indicator  of  physiological  state  (Bell,  1985;  Langdon,  1987).  Features  of 
this  instrument  include:  1)  a  quartz  culture  chamber  extending  the 
response  spectrum  into  the  ultraviolet;  2)  a  system  to  control  the  pH  (total 
C02)  and  dissolved  oxygen  that  allows  unattended  deployments  in  situ  up 
to  48  hours;  3)  a  nicroprocessor  based  control  and  data  acquisition  system 
providing  flexibility  in  sampling  through  software  modification. 

An  Ocean  Optics  spectrometer,  purchased  with  funds  from  this  ONR- 
SORP  grant,  was  added  to  SUPA.  The  present  configuration  uses  a  dual 
fiber  optic  unit  with  two  SUPA  housings.  One  fiber  is  incorporated  into  the 
SUPA  housing  that  contains  a  culture  of  the  phytoplankton  species  of 
interest;  the  second  fiber  is  incorporated  into  an  identical  SUPA  housing 
that  contains  DI  wa:er.  The  readings  from  the  two  fibers  simulate  a  dual 
beam  spectrophoton  eter  and  provide  information  on  the  spectral 
absorbance  of  the  phytoplankton  culture  in  real  time  as  the  SUPA  units  are 
moved  throught  the  upper  mixed  layer.  This  capability  used  with  the  other 
SUPA  probes  poten Tally  provides  the  ability  to  monitor  the  short  time 
scale  changes  in  quantum  efficiency  exhibited  by  the  enclosed 
phytoplankton  cells  that  have  a  known  PPFD  history  in  a  natural  water 
column. 


MODELING 

Kamykowski  (1990),  using  a  random  walk  simulation,  investigated  how 
physical  and  biolog  cal  vectors  combine  in  a  three-stratum  mixed  layer  to 
influence  phytoplan  cton  distributions.  The  average  vertical  eddy 
diffusivity  coefficient  scaled  after  Sundby  (1983)  decreased  by  one  order 
of  magnitude  in  ea(  h  successive  stratum.  Following  the  discussion  in 
Denman  and  Gargeit  (1983),  the  upper  layer  included  Langmuir  circulation 
as  a  large  turbuleni  eddy.  Two  different  motility  patterns  based  on 
buoyancy  (sinking,  neutral  and  floating  cells)  and  swimming 
(dinoflagellate,  slow  ciliate  and  fast  ciliate)  were  compared  under  various 
physical  forcings  ('^  amazaki  and  Kamykowski;  1990).  They  also  examined 
population  distributions  of  variously  motile  phytoplankton  but  in  a 
continuous,  unstrati 'ied,  wind-driven  mixed  layer  characterized  by  an 
exponentially  decreasing,  dissipation  profile  with  a  lower  boundary  defined 
by  the  Ekman  depth.  The  PPFD  exposure  of  each  cell  was  recorded. 
Distributions  of  a  c  inoflagellate  at  three  different  wind  speeds  and  the 
profiles  of  the  PPFD  exposure  for  each  cell  differed  markedly  from  that  of 
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a  cell  in  the  absence  of  turbulence  and  to  a  threshold  value  below  which 
photo-stress  is  assumed  to  occur. 

The  modeling  jffort  during  the  term  of  the  ONR-SORP  grant  emphasized 
incorporation  of  a  photoadaptaton  submodel  into  the  previously  described 
random  walk  Lagraiigian  models  (Kamykowski,  et  al.,  1994)  and  into  a 
version  of  the  Price  et  al.  (1986)  Eulerian  model  of  the  upper  mixed  layer 
(Janciwitz  and  Kamykowski,  1991).  A  model  by  Denman  and  Maira  (1986) 
provides  a  time-dej  endent  response  to  the  PPFD  history  experienced  by  a 
phytoplankton  cell.  The  completed  Lagrangian  turbulence-behavior- 
photoadaptation  model  provides  the  depth  and  instantaneous 
photosynthetic  rate  of  an  individual  algal  cell  with  an  assigned  sinking  rate 
that  is  randomly  mixed  based  on  an  assigned  wind  speed  in  an 
independently  randon  PPFD  field.  Successive  runs  of  the  model  provide  a 
population  ensemble  of  the  species  under  consideration  that  reflects  the 
range  of  photosynth.‘,iic  history  of  a  group  of  cells  at  a  given  time  and 
place  in  the  water  column.  Alternate  computation  routines  are  under 
investigation  (Yamazaki,  personal  communication).  The  Eulerian 
turbulence-convectic  n-photoadaptation  model  calculates  the  changes  in 
the  temperature  anc  the  photosynthetic  rate  profiles  over  the  diel  cycle. 
The  phoiosyntheiic  rate  at  a  given  depth  is  influenced  by  the 
redistribution  of  the  phytoplankton  cells  due  to  vertical  mixing  and  to  the 
average  effect  of  eisemble  photoadaptation  of  the  mixed  population. 
Comparisons  with  published  field  data  are  reported  (Kamykowski  et  al,, 

1996) .  These  workiig  models  and  their  ongoing  refinements  provide  a 
framework  and  motivation  for  the  specific  field  experiments  planned  in 
the  future.  Comprel  ensive  data  sets  are  required  to  fully  test  model 
parameterizaiions  <nd  performance. 

Additional  mo  ieling  has  examined  dinoflagellate  orientation  in  a 
Lagrangian  context  ^Kamykowski,  and  Yamazaki,  1996;  Kamykowski  et  al., 

1997) .  The  backgromd  for  these  models  is  summarized  in  Jones  (1993) 
and  Kamykowski  (1995).  In  general,  the  the  individual  cells  adjust  their 
orientation  based  0)i  their  environmental  experiences  (ex.  PPFD,  nutrients, 
temperature)  in  the  upper  ocean  and  on  their  internal  physiological  state 
(ex.  storage  pool  si:.e,  cell  cycle  stage).  These  models  provide  a  different 
and  possibly  more  ealistic  result  compared  to  cells  that  simply  swim 
toward  the  light  du  ing  the  daytime  and  away  from  the  surface  at  night. 
Supporting  laboratory  experiments  are  underway 
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DISCUSSION 

Although  add  tional  funds  were  not  requested  through  the  ONR-SORP 
program  after  199^  due  to  the  shifting  research  interests  of  Dr.  Tom 
<^urtin,  the  work  described  above,  to  which  the  ONR-SORP  funds 
conmbuted,  is  con  inuing.  SUPA  deployments  in  support  of  modeling 

destaWe 
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